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Tjlj- , Abstract. A generalization of determinant formulas for the classical solutions of 

(N ■ Painleve XXXIV and Painleve II equations are constructed using the technique of 

Darboux transformation and Hirota's bilinear formalism. It is shown that the solutions 
^ . admit determinant formulas even for the transcendental case. 
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. 1. Introduction 

It is well-known that solutions for the Painleve equations play a role of special functions 
in the nonlinear science . Originally, Painleve derived these equations in order to find 
new transcendental functions determined by second order ordinary differential equations 
possessing so-called the Painleve property. Recently, irreducibility of solutions of the 
Painleve equations has been proved by Umemura et al §|. Umemura first gave a 



O 

C/2 



H ' rigorous definition of classical functions: starting from the field of rational functions, 
if a function is obtained by finite numbers of iterations of the following permissible 
operations, 

• differentiation, 

• arithmetic calculations, 

• solving homogeneous linear ordinary differential equation, 

• substitution into Abelian functions, 

then that function is called "classical". 

Umemura proved that solutions of the Painleve equations are not classical in general 
in the above sense. However, it is known that they admit classical solutions for special 
values of parameters. One is the rational or algebraic solutions, and another is the 
transcendental classical solutions which are expressed by rational functions in various 
special functions and their derivatives. 
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In this article, we discuss the Painleve II equation (Pn), 

^ = 2„'-4.. + 4(a. + i). (1) 
and the equation of No. 34 in Gambier's classification |Q, 

which we call the Painleve XXXIV equation (Pxxxiv)- In equations (|lD and (0), we 
adopted diff'erent scale from their canonical forms 

— = 2u^ + su + (3, (3) 

~ ( rfj ) ~ ^ +1 = 0, (4) 
for the clarity of expression of solutions. The scale transformations. 



2. 

lead the canonical forms to our ones. We denote Pii[a] and Pxxxiv 

respectively, when it is necessary to show the value of parameter explicitly. These 
equations are related to each other by the Miura transformation |^, 

w = — --u^ + 2z , (6) 

dz 

and its complement, 

dw 



U 



dz (7) 



2w 

Umemura et al also proved that Pn (and thus Pxxxiv) admits transcendental classical 
solutions which are expressed by the Airy function when a is an integer, rational 
solutions when a is a half-integer, and that otherwise the solutions are non-classical. 

It is known that the Painleve equations (except for Pi) admit the Backlund 
transformations(BT) which form the affine Weyl group ^, |TT|. For example, 
BT of Pn is given by, 

Aot 

S : S{u) = u + , 5(0;) = —a, (8) 

-2z 



dz 
~du 



T : T{u) = -u+ , ^ , T{a) = a + 1, (9) 



— -u'^ + 2z 
dz 



and < S,T > forms the affine Weyl group of type A^i \ Starting from a suitable "seed" 
solution, we obtain "higher" solutions by applying BT to it, which are expressed by 
rational functions in the seed solution and its derivatives. 

It is interesting and important to note two points on classical solutions for the 
Painleve equations. The first is that, in general, the classical solutions are located 
on special points in the parameter space from the viewpoint of symmetry. Namely, 
transcendental classical solutions are on the walls of the Weyl chambers, and rational 
solutions on their barycenters. The second is that the classical solutions have additional 
structure. Namely, it is known that some of the classical solutions admit such 
determinant formulas that they are expressed by log derivative of the ratio of some 
determinants. Moreover, application of BT corresponds to increment of size of the 
determinant. In fact, both Airy function type and rational solutions for Pn and Pxxxiv 
admit such determinant structure. 

Then, does such determinant structure exist even for non-classical solutions, or are 
non-classical solutions so "transcendental" that do not admit even such structure? Here 
we also note that such determinant structure are shown to be universal among the 
soliton equations due to the celebrated Sato theory |13|. 



The purpose of this article is to generalize the determinant formulas of the classical 
solutions of Pii and Pxxxiv, and show that such determinant structure is universal 
among the solutions of Pn and Pxxxiv- The key to this generahzation is to use the 
technique of the Darboux transformation together with Hirota's bilinear formalism. 

This article is organized as follows. In section ^ we summarize the derivation of 
Pxxxiv and classical solutions of Pn and Pxxxiv- We illustrate our basic ideas in section 
^. Our main results are presented in section ^. The proof for our results are given in 
section ^. Section |^ is devoted to summary and discussions. 

2. Pxxxiv and classical solutions of Pn and Pxxxiv 

It is well-known that Pn is derived as the similarity reduction from the modified KdV 
equation 0, 

Ut + ^U^Ux - ^Uxxx = 0. (10) 

In fact, putting 

equation (|10|) is reduced to equation ([1|). Similarly, Pxxxiv (@) is derived as the similarity 
reduction from the KdV equation, 
3 1 

wt - -ww^ - -w^xx = 0, (12) 



by putting 



w{x, t) 



(w(z) -2z) 



X 



It is known that Pn (|T]) is bilinearized to 

-f,)g.f = 0, 

Dl + {Az - 3iJ,)D, - 4 (^a + 

by the dependent variable transformation, 
d , g 

'' = d-z^''^r 



(3t)i/3- 



g-f = o, 



(13) 

(14) 
(15) 

(16) 



where fi is an arbitrary function in z, and Dz is Hirota's bihnear differential operator 
defined by 

d d ^ ^ 



D^g-f 



dz dz' 



gi^)fi^'] 



(17) 



The bilinear equations ([liD and ([T5| ) are regarded as those for Pn, but it is also possible 
to derive Pxxxiv as follows. First we divide equations ([l^ ) and (|15|) by p. Then, 
putting 

g_ 

r 



w 



2— log/-/i + 2z. 



and using the formulas |16, 17 



D.g- 


f 


P 




Dig- 


f 


P 




Dig- 


f 


P 



dz ' 
fd^ 
[dz^ 



+ w + fi-2z]il), 



we obtain 



o/ ^ ^ d 

— + 3iw + fi-2z) — 
dz-^ dz 



' d dw , \ , 

2w- 4a V = 0, 

dz dz / 



dz"" 



+ w- 2z\il) = Q. 



(19a) 
(19&) 
(19 c) 

(20) 

(21) 



Compatibility condition of the linear equations (pOf ) and (^) gives Pxxxiv 

Let us discuss the determinant expressions for the classical solutions of Pn and 
Pxxxiv [|10], |13|, |T5[ . In the following, the determinant with size zero should be regarded 
as 1. 



If we choose /i = then it is known that the bihnear equations (|T^) and (|T5|) are 
satisfied by 

g = PN^x. f = pN, a = N, NeZ>o, (22) 



Pn 



fiO) fil) 
/(I) /(2) 



jiN-1) j!{N) . . . j{2N-2) 



Po = 1, 



(23) 



where /("^ 



dz 



f and / satisfies the Airy equation, 
,=2zf. 



Thus, the Airy function type solutions of Pii[A^] and Pxxxiv[^] are given by 

d' , Pn+i 
u = -^hg , 

dz 

d"^ 

w = 2— log Pn, 



(24) 

(25a) 
(256) 



respectively. 

The rational solutions are obtained for the case of /i = 0, and there are two 
expressions for them. One is the Schur function type expression given by 



g = (Ttv+i, f = (Tn, a = N + -, N e Z>o, 



Qn Qn+i ■ ■ ■ Q2N~1 
<iN-2 QN-I ■ ■ ■ (l2N~3 



(26) 



(27) 



Q-N+2 Q-N+3 ■ ■ ■ Ql 

where are the polynomials in z defined by 

E Ikiz) = exp (z^ + \) , Qkiz) = for A; < 
Another expression is given by Hankel determinant as 

g = kn+1, f = i^n, a = N+^, NeZ>Q, 



Km 



ao ai 
ai a2 

(IN—1 CiN 



ajv_i 

a-2N-2 



Ko = 1, 



(2^ 



(29) 



(30) 
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where a„, n = 0, 1, 2, ■ ■ ■ are the polynomials in which are recursively defined by 
da 1 

an = "~ + V akan-k-2, n > 0, ao = z. (31) 

fc=0 

Thus, the rational solutions for Pii[A^ + |] and Pxxxiv[^ + ^] are given by 

u = —\og = i-log > (32a) 

az (Tat az kn 

£ (f 
w = 2z + 2—\ogaN = 2z + 2—\ogKN, (326) 

respectively. 

We note that it is possible to generalize the above result for negative a by the 
refiection symmetry of Pn and Pxxxiv, 

u{—a — 1) = —u{a), (33) 

and 

w{—a) = w{a), (34) 
respectively. Indeed, this symmetry is generated by ST in equations (|^) and (^. 

3. Darboux Transformation 

The technique of Darboux transformation is well developed in the soliton theory For 
example, in the case of KdV equation (0), we start with its auxiliary linear problem. 



3 ^ 3 ^ (35) 



where A is the spectral parameter. Then, one can show that equations (|35D are covariant 
with respect to the Darboux transformation \E' — > \E'[A^], w w[N] defined by 

' ' lVr(*i,f2,...,f«) ^ ' 

w[N] =w + 2— logT^r(^i, ^2, ■ ■ ■ , ^iv), (37) 

where is a solution of the linear equations (^) with A = A^, and Wr is a Wronskian 
with respect to the indicated functions. Thus, choosing an appropriate seed solution of 
KdV equation as w, one can construct series of exact solutions by using this method. 
For example, starting from the solution w = 0, then we obtain the Wronskian expression 
of A^-soliton solution Although most of the expression would be somewhat formal 
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except for the soliton and rational type solutions or their variants, it is important that 
we can express quite wide class of solutions in terms of determinant. 
The above result is recovered by the following bilinear equations: 

D^Dt - ^-Dl - \wd)j F-F = 0, (38) 

[dI + X + w)G-F = 0, (39) 
Dl - 4Dt + 3(-A + w)D^] G ■ F = 0, (40) 

where F = l^r(^i, ^'2, ■ ■ ■ , ^Af), G = Wr{'^i,'^2, ■ ■ ■ n,"^)- We omit the detail, 

but we can directly prove that these bilinear equations holds as the identity of the 

determinants from the linear equations (|35|). Indeed, we recover "usual" bilinear 
equations by putting w = 0, 

D,Dt - ^Dt^ F-F = 0, (41) 

DI + \)G-F = 0, (42) 
Dl - ADt - 3XD^ G-F = 0. (43) 



^From these bilinear equations, we recover the KdV equation ([1^ and the auxiliary 
linear problem (^) by putting as 

«^ = 2^1og^, * = F ^^^^ 
Now, keeping this correspondence between the Darboux transformation and the 
bilinear formalism in mind, let us illustrate our strategy to generalize the determinant 
formulas for the solutions of Pn and Pxxxiv- We have the bilinear equations ( |T^ ) and 
(P^) for Pxxxiv- Clearly, these bihnear equations correspond to equations (^) and 
(^). Thus, if we could obtain "generalized" bilinear equations which correspond to 
equations ( ^9l) and (^0|) , it might be possible to find the Darboux transformation for 
Pxxxiv and Pn which leaves the linear equations (|20|) and (|2l|) covariant. Then we 
obtain a generalization of the determinant formulas. 
In the next section, we present our main results. 

4. Main Results 

In this section, we present a generalization of determinant formulas for the solutions of 
Pii and Pxxxiv- 

We start with a pair of linear equations (|20| ) and (pT]). As we mentioned in 



section these equations are compatible provided that w satisfies Pxxxiv [«]- In other 



words, equation (|21| ) follows if w and ip are solutions of Pxxxiv[ci;] and equation (|20|), 
respectively. Then we have the following theorem. 

Theorem 4.1 Let w he a solution of i^xxxivlt^]; one? ip^ he a solution of the linear 
equation i\2di). We define two sequences ipn, = 0, 1, 2, ■ ■ ■) hy 

w 



W 



dz 2ipQ 

dipn-l 



n-2 



dz 



i>ki)n-2-k, n > 0, 

fc=0 
n-2 

+ V'O X! fk'Pn-2-k, n> 0. 



(45) 
(46) 
(47) 



fc=0 



We define Hankel determinant , N & Z hy 

det(V'j+j-2)ij=i,-7v > 
tn= <^ 1 = . 

N <0 



Then, 



satisfies 



N 



2W 



d dW 



dz dz 



4(a + N) 
0, 



N 



where 



W = W + 2— log Tat. 



(48) 

(49) 

(50) 
(51) 

(52) 



Thus, W satisfies Pxxxiv[a + N]. 



Similar formula for Pn is obtained by applying the Miura transformation (P). The linear 
equations (pO|) and (|2l|) are reduced to 



Then we have the following theorem. 



(53) 



Theorem 4.2 Let u he a solution of Pii[q\, and r^r he Hankel determinant given hy 
equation ^^), where ipn and ipn are defined recursively hy 

du 9 



Then, 



'4'v 



U 



dtp, 



n-l 



dz 

dip n-l 



n-2 
k=0 



n > 0, 



n-2 



dz 



k=0 



d tn+i 
— log 

dz tn 
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(55) 
(56) 

(57) 



satisfies Pii[a + A^]. 

As mentioned in the previous section, Theorems O and if! is the direct consequence 
of the following proposition: 



Proposition 4.3 The following bilinear equations hold. 
{D1 + w - 2z)rN+i ■ rjv = 0, 
Dl + {3w - 2z)D, -A(a + N +\ 



TN+I ■ Tat = 0. 

In fact, dividing equations (|58D and (^) by and using the formula [|16], 

DzTn+1 ■ Tn _ d^^ 

dz ' 



N 



T. 



N 



dz'^ 

d^ 



N, 



(58) 
(59) 

(60a) 
(606) 
(60 c) 



we obtain the linear equations (|50|) and (0) in Theorem [4.1| . Similarly, dividing 
equations (^8]) and (|59| ) by t^+iTn and using the formulas |jl6|, p!7| , 

DzTn+1 ■ Tn 



tn+iTn 



DItn+1 ■ Tn j^jj2 

tn+iTn 
D'Itn+i ■ Tn d'^U 



with 



we get 



Tn+iTn 



d^ 



dz^ 



+ WV + W 



V = ^ log {tn+i'Tn) 



V + U^ + w-2z = 
d'^U 



dz^ 



+ 3UV + {3w - 2z)U - 4 a + AT + 



(61a) 
(616) 
(61c) 

(62) 

(63) 
(64) 
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Eliminating V, we obtain 

^ = 2U^-4:zU + 4:(a + N + -]. (65) 
dz^ \ 2 J 

These results imply that even non-classical solutions possess the same determinant 
structure as the classical solutions. Indeed, the known determinant expressions for the 
classical solutions are recovered as special cases. In fact, starting with a solution w = 2z 
for Pxxxiv[|]) the linear equation ( 120|) yields 

(4 - * = °' 

from which we get ipQ = z without losing generality. Then, we have from the recursion 
relations (H), (|6|) and (0), 

= —1 h 2^ llJk^pn-2~k, n> 0, ll)Q = Z, (67) 

'^^ fc=0 



- + z^ ipk^Pn-2-k, n > 0, v^o = 1- (6^ 



dz 

fc=0 



Thus, we have a series of rational solutions given by equations (^8]), (|67D, ( |68D and 

^ = ;flog^, (69a) 
d^ 

W = 2z + 2—\ogTN. (696) 

The case of > agrees with the Hankel determinant expression for the rational 
solutions discussed in section 0. 

Next, noticing that w = is a solution of Pxxxiv[0], we find from equation (pl|) 
that ipQ is determined by 



^-2z)^ = 0, (70) 



and the recursion relation (^) with equation ( ^5]) is reduced to 

= — 1 • (71) 

dz 

Thus, we have a series of solutions given by equations (|48|), (|70|), (|7lD and 

f/ = ;flog^, (72a) 
dz tn 

W = 2-^\ogTN, {72b) 
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for > 0, which agrees with the Wronskian expression for the Airy function type 
solutions discussed in section 0. However, the Airy function type solutions are "singular" 
case in our formula. From equations ( |45| ) and ( ^Tj) , we find that (pn = for all n and thus 
Tiv = for < 0, which does not give meaningful result. We note that this phenomena 
is related to the symmetry of the r sequence which will be mentioned in section ^ and 
the correct result for iV < case is obtained by virtue of symmetry. We also note that 
this is the only singular case and our determinant formula works for other cases. 



5. Proof of theorems 



In this section, we give the proof of the Theorems f4.1| and |4.2| . Since we see that these 
theorems follows immediately from Proposition it is sufficient to prove it. We have 
to prove both > and < cases, but we concentrate on the former case, since the 
latter case is proved in similar manner. 

The bilinear equations (^) and (|59| ) are reduced to the Pliicker relations, which 
are quadratic identities of the determinants whose columns are shifted. Thus, we first 
construct such differential formulas that shifted determinants are expressed by operating 
some differential operator on the original determinant. For this purpose, we introduce 
the following notation: 



Definition 5.1 Let Y he a Young diagram Y 
N X N determinant t^y by 

V'o V^i ■ ■ ■ i'N-h-l i'N-h+ih 
^1 ■■■ -ipN-h tpN^h+l+ih 



tny 



WN ■ ■ ■ '^2N-h-2 W2N-h-l+ih 

Then, we have the following differential formulas 



^2 



(ii, ^2, ■ ■ ■ ; "^/i)- Then we define an 



'^N-2+i2 V'TV-l+ii 



2Af-3+i2 4'2N-2+i 



(73) 



Proposition 5.2 

d 

dz 

( d^ w\ 



w 



+ 2Nz tn 



T Nrm — Tm 



( d^ 3d ldw\ 

w \ d , 

-- + 2Nz)— + 2{a-l + N) 
2 J dz 



(74) 
(75) 
(76) 
(77) 

(78) 
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Tat^ — TiV|ip + TnI 



d 



1 dw 



2z— + 2N^ + AN{a - 1) - --j- - 3(a - 1 



dz 



4 rfj; 



(79) 



The proof for Proposition p.2| is an important step. However, since this requires 
straightforward but tedious calculations, we will give it in [Appendix A . 
Finally, we prove Proposition [4.3| . /^From the Pliicker relations we have, 

TAT+igTAr - TTv+ipTjvp + TAr+irjv^ = 0, (80) 
'TN+1\P'Tn — '^N+lm'TNu + TN+lTNvm = 0. (81) 

]) and (p3) from equations (BH) 



By using Proposition |5.2| , we get the bilinear relations 
and (|81|), respectively, which is the desired result. 



6. Summary and Discussions 

In this article, we have presented determinant formulas for the solutions of Pxxxiv 
and Pn which are valid also for non-classical solutions by using the technique of the 
Darboux transformation and bilinear formalism. The solutions of Pxxxiv + A^] and 
Pn[Q; + A^], N E Z, are expressed by determinants whose entries are constructed from 
the solution of some linear equations. Moreover, coefficients of those linear equations 
include the solution of Pxxxiv [«] and Pii[a], respectively. We have also shown that 
known determinant expressions for classical solutions are recovered as special cases. 
This result implies that determinant structure of the classical solutions is universal 
among the solutions of Pn and Pxxxiv- 

Finally, let us discuss the relation with the Toda equation, which will be a key for 
generalization to other Painleve equations. 



In general, the r function for Pn is introduced through its Hamiltonian [10 

1 



by 



Hii{v^u,z]a) = -^v + {—u + 2z)v — Aau, 



Hu{v,u,z;a) = -^\ogT{a). (82) 



We note that we obtain Pn [a] for u from the canonical equation, 

du dHii dv dHii 
dz dv ' dz du 
Then it can be shown that m(q;), which is a solution of Pii[q;], is expressed as 



^3) 



d , ria + 1) , , 

u{a) = -\og \ \ 84 
dz T{a) 
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By applying BT, we obtain a sequence of r functions {Tj\f}]\f^z, where = r(a + A^). 
Okamoto has shown that BT of Pn is governed by the Toda equation on the level of r 
function: 



Proposition 6.1 (Okamoto /(7^/^ satisfies the Toda equation, 

— log Tat = ctv 2 ' 

where cn is a non-zero constant. 

It is well-known that the solution of the Toda equation ( pS] ) is expressed by 



^^ = det(^-^/ , (86) 

i,j=l,---N 



where / is an arbitrary function in z. Equation (|86|) is sometimes referred as "Darboux's 
formula". However, one should note that Darboux's formula (|86|) is valid under the 
condition, 

ro = l, ri = /, N>0 (87) 



(In practice, tq can be a constant.) In [jlO|, it is pointed out that the Wronskian 
expression for the Airy function type solutions of Pn is a consequence of the Darboux's 
formula. We demonstrate how we could apply the Darboux's formula for this case. We 
have a solution of Pn for a = 0, 

M(0) = ^logV^, ^^ = 2z^. (88) 
Therefore, we can choose the r functions as 

r(0) = To = 1, r(l) = n = ^. (89) 
Then we have a r sequence for the Airy function type solutions, 

• • ■ r_3, r_2, r_i, Tq = r(0) = 1, Ti = ip, Ts, rg • • • . (90) 

Now, as was mentioned in section |^, we have a reflection symmetry (|33|) on u which 
implies a symmetry on the r function, 

r(-a) = r(a). (91) 

In the case of the Airy function type solutions, this symmetry induces a symmetry on 
the r sequence as 

r-N = TN. (92) 
By virtue of this symmetry, we see that the r sequence is divided into two parts as 

■ ■ ■ ^3, ^2, Ti, To = r(0) = 1, Ti = i), T2, T3 ■ ■ ■ . (93) 
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Fortunately enough, since it can be shown that this r sequence is governed by the Toda 
equation (|85D with = 1, we could apply Darboux's formula for > and < 
separately. We note that if we prolong the r sequence for < following to the Toda 
equation (p5| ) with the condition { ^7\ ) without taking the symmetry into account, we 
have Tat = for < 0. This corresponds to the "singular" phenomena mentioned in 
section ^. 

In the case of rational solutions, the r sequence is again separated into two parts. 
However, we cannot apply the Darboux's formula to this case. Let us take a solution of 
Pii for a = 1/2, 

We can choose the r function as 

t(1-]=to = 1, T(^]=n = z, (95) 



.2/ V2. 

and we have a r sequence for the rational solutions. Now the symmetry (0) implies a 
symmetry on the r sequence, 

tn = r_iv_i, (96) 

and thus we have, 

■■■^2, Ti, To, To=t(]-)=1, Ti = Z, T2, Tg ■ ■ ■ , (97) 



.2 

and again the r sequence is separated into two parts. However in this case, these r 
functions are shown to satisfy the Toda equation of the form, 

^logr^ = — ^-z. (98) 

Thus, we cannot apply Darboux's formula for the rational solutions. Of course, by 
introducing a gauge on the r function as 

(TN = e^'/' TN, (99) 
then ajy satisfy the Toda equation 

^loga. = ^^. (100) 

However, now ctq is not a constant. 

There are two points for being able to apply the Darboux's formula. The first is 
that T sequence should have a symmetry which is induced from the Painleve equation 
itself so that the r sequence is separated into two parts, which is necessary to apply the 
Darboux's formula without inconsistency. The second is that r sequence should satisfy 
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the Toda equation of the form ( pSf ) under the condition (|^). Both conditions are 
satisfied for the Airy function type solution , but the second condition does not hold for 
rational solutions. Now, for solutions which corresponds to generic a, the symmetry on 
the T function ( [91| ) does not induce any symmetry on the r sequence. This observation 
shows that it is not trivial that generic r function admit determinant formula, even if 
it satisfies the Toda equation. 

Despite of unavailability of the Darboux's formula, we could obtain the determinant 
formula for rational solutions. This is due to the determinant formula of the general 
solution of Toda equation of C-type (Toda equation with the symmetry (p6l)) obtained 
in [H. 

Conversely, the general determinant formula for Pn strongly implies that it is 
possible to construct the determinant formula for the general solution of Toda equation 
in general setting, i.e., the Toda equation admits the determinant formula for the general 
solution without any symmetries or gauge on r sequence. Then, our results might be 
regarded as the special case of such general solution for the Toda equation. 

Moreover, it is known that the Backlund transformations for the Painleve equations 



(except for Pj) is governed by various types of Toda equations p|, |^, [1^, |n[]. Thus, 
it might be possible to present general determinant formula also for other Painleve 
equations. We shall work out this problem in the next publication. 
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Appendix A. 



Here, we give the proof for Proposition 5.2 



We first prove equation (fflD . Notice that ttvd is expressed by 



TNd 



( ^1 

^^2 



^2 



Af+1 



/ A 



11 



A 



A21 A 



12 



22 



A 



N2 



Aatat / 



(Al) 



where Aij is the (z, j)-cofactor of tn and A ■ B denotes a standard scalar product for 
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N X N matrices A = (aij) and B = (bij) which is defined as 



N 



A - B = aijbij = trace A ^B. 



(A2) 



The first matrix of equation ([Al| ) is rewritten by using the recursion relation (^61) as 





If, 


dz, 
d 

dz 


dz 


dz 


d ' 

dz 


d 

dz 


d ' 

dz 



w 



2^0 



^0^1 + ^1^0 



N-2 N-1 
k=0 k=0 



The above matrix in the second term is separated as 











N-2 \ 
i)k^N-2-k 

k=0 
N-1 

J2 i'ki'N-i-k 

k=0 

2N-3 

Y ^k'ip2N-3-k 
k=0 



\ 

1pN-l1pO 



+ 



N-2 N-1 

Y i^ki^N-2-k i^kipN-i-k ■ ■ 

k=0 k=l 

I N-2 ^^ 

^0 ■ ■ ■ I] iJki'N-2-k 

k=0 
N-2 

^0^1 ■ ■ ■ Y i^ki'N-i-k 

k=0 



2N-3 

Y ^k^2N-3-k 
k=N-l 



tpo^N- 



N-2 

Y ^k1p2N-3-k 
k=0 



(A3) 



(A4) 



Each of these terms gives zero contribution in equation ( |A1| ). Hence we obtain equation 
(0)- 



Next we prove equation (^). We consider 
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/ ^1 ^"2 



( A, 



oil Api2 
^□21 Aq22 



•□Afl Apjv2 



■ ApjvAT y 



(A5) 



where Apjj is (i, j)-cofactor of r^vp. The first matrix in the right-hand side is equal to 



d 


d 

If, 


d 

■ -rVN-2 
If, 


d 


dz 


dz 


dz 


dz 


d ' 

dz 


d 

dz 


d ' 

■ -rW2N-3 

dz 


d 

-rW2N-i 

dz 



w 



^0 



N-2 

1pk1pN-2-k 
\ k=0 



W 



2^0 



ro 

V^oV'i 

tpoi^N- 





N-l 

k=l 
N-3 

J2 i^ki^N-s-k 

k=0 
N-3 

i)kijN-2-k 

k=0 
N-3 









^N'ijjQ 



2N-A 2N-2 

Y '^ki'2N-4-k Y '^ki'2N-2-k 

k=N-2 k=N 

N-l 



Y ^k1pN-l-k 

k=0 

N-l 

Y i'ki^N-k 

k=0 



N-l 



Y ^k'ip2N-4-k Y '^k^2N-2-k 
k=0 k=0 



(A6) 



Taking the scalar product, the first and second terms give — T^rp and — r^r, respectively, 

dz 2 

and the third term vanishes. Hence we have equation (|7S|). 

Next we prove equation (^). We consider the following equality. 



/ i^2 

^3 



^4 



i^N+1 \ 
i^N+2 



\ IpN+l 'ipN+2 



( A 



11 



A 



A21 A 



12 



22 



%l)2N j V Aati AAr2 



Aijv \ 
A2JV 

Aatat ) 



(A7) 



The first matrix of riglit-liand side is rewritten as 



d 

-1-^1 

If, 

dz 


d 

-1-^2 ■ 

If, 

dz 


d 

dz 


d 

dz 


dz 


d ' 

■ -r'ip2N-i 
dz 



w 



2^0 











N 



N-1 

V k=l k=2 



W 



2^0 



^0^1 + ^iV'o 

^0^2 + V'lV'l 



V 



Here we note tliat ipnS also satisfy 

dipn 







2N-2 

Y ^fc^2Af-2-fc 
k=N 

N-1 \ 

J2 i^kipN~i-k 

k=0 

Y i^k'ipN-k 

k=0 



N-1 

Y ^k^2N-2-k 
k=0 



dz 



dw 

— 4a n-2 

{2Z - W)i)n-1 + 2{n - l)i)n-2 - Y] 1pk^n-2-k, 



wliicli is proved by induction from equations (pOD , (l2l| ) and (|46|) . The first term 
right hand side of equation ( |A^ ) is rewritten using equation (|A9|) as 



(2^ - w) 



/ V'o i^i 

^1 ^2 



'ipN 



+ 



\ i^N-l i^N ■■■ i^2N-2 J 



2^0 2^/^! 

V 2(iV-l)^^_2 2(iV - l)V'7v-] 



\ 

2^N-1 
2{N - l)^27V-3 J 
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+ 



/ 2^0 
2^1 



2(iV-l)^^_2 \ 
2(iV-l)V';v-i 



V 2^N-l ■■■ 2{N- l)^2N-3 J 

/ 

^0 ^1^0 



dw 

4a 

dz 



N-2 

J2 ^ktpN-2-k 
\ k=0 



( 



dw 

dz 



- Aa 



4V^o 



^2 



IpQlpi 



i^QipN- 



N-l 

^k^N-i-k ■ 

k=l 
N-2 

J2 lpk'lpN-2-k 

k=0 

N-2 

J2 i^ki'N-l-k 
k=0 

N-2 

Y '^k'ip2N-3-k 
k=0 







2N-3 

Y ^ki'2N-3-k 
k=N-l 



(AlO) 



Applying the scalar product on these terms, we obtain equation (^). We get equations 
([r7|)- ([79|) by similar calculations. 
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